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Spiral waves are a ubiquitous feature of the nonequilibrium dynamics of a great variety of 
excitable systems. In the limit of a large separation in timescale between fast excitation and slow 
recovery, one can reduce the spiral problem to one involving the motion of a free surface separating 
the excited and quiescent phases. In this work, we study the free surface problem in the limit of 
small diffusivity for the slow field variable. Specifically, we show that a previously found spiral 
solution in the diffusionless limit can be extended to finite diffusivity, without significant alteration. 
This extension involves the creation of a variety of boundary layers which cure all the undesirable 
singularities of the aforementioned solution. The implications of our results for the study of spiral 
stability are briefly discussed. 


Understanding the behavior of spiral waves in excitable media remains an important and challenging problem [Q. 
From the perspective of numerical simulation, two component reaction-diffusion models have been shown to 
capture the important features of spiral patterns, in particular the transition from rigid rotation to meandering and 
the phenomenology of the nonlinear meandering state QJ|] . More recently, simulations Q , and a numerical stability 
analysis Jij] based on the free boundary limit [|| (see below) have confirmed that the finite thickness of the front 
(separating the excited region from the quiescent one) is not crucial for the meandering behavior. It is therefore of 
interest to pursue analytical techniques (which make use of this free boundary limit) in the hope of gaining a more 
fundamental understanding of the nature of this instability. 

This work reports progress towards the aforementioned goal of having an analytic theory of spiral waves. Specifically, 
we revisit an approach due to Pelce and Sun who derived a spiral solution in the case of zero diffusion for the 
“controller” variable. Their solution exhibits singular behavior near the spiral “tip” (for example, the front curvature 
has a discontinuous derivative across the tip point) which raises questions regarding the validity of the solution and 
to date has precluded a full stability analysis. Here, we show how the inclusion of small but finite diffusivity provides, 
via the introduction of boundary layers, for a regularization of the singular behaviors. This therefore confirms the 
physical validity of their construction. The implications of our results for a (future) calculation of spiral stability are 
discussed at the end. 

The free boundary approach starts from the equations coupling a concentration field v(x, t) to an interface between 
an excited state of the medium (“+”) and the quiescent state (“-”). In so-called “Fife”-scaled units |ll|| , 

=9±~ RV± + DN7 2 v (1) 

where /i is a positive constant, g+, c/_ are positive, negative constants and v± refers to the field in the +,- regions 
respectively. The field obeys the boundary conditions at the interface 

Cn T ^ — Vint (2) 

where c n is the normal velocity, n the curvature, and the value of v (as well as its normal derivative) is continuous 
across the interface. 

This system of equations is not rigorously derivable from the original reaction diffusion system 


du ~ f(u,v ) 

at e 

(3) 

d'v 

— = D\7 2 v + g(u,v) 

(4) 
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The reason for this is that replacing the “propagator” field equation for u(x, t ) by an interfacial boundary condition 
is only valid asymptotically as e —> 0. The coefficient g. of the linear term in eq. (1) is formally small, of order e 1//3 . 
Hence we cannot rigorously keep this linear term without keeping additional terms of the same order as well. As we 
have already mentioned, however, simulations show that eqns. (1) and (2) do capture the spiral phenomenology of 
interest, at least for finite diffusivity D. Furthermore, an exact numerical steady-state solution 0 and subsequent 
stability analysis ji]] directly supports this conclusion. On the other hand, results obtained by dropping the linear 
term and working at small D seem at present to be rather unphysical |u|[7j]. We will therefore adopt eqns. (1, 2) as 
our fundamental model and proceed to consider the small D limit. 



FIG. 1. A typical spiral, with rt, 0/(r), 9b(r) indicated 


This free 
field v(x, t) 


boundary problem 
—> v(r, 9 — ujt) one 


with D = 0 was first tackled by Pelce and Sun . If one assumes a uniformly rotating 
obtains 


dv± 

-“-Q0=9±~» v 


( 5 ) 


This can be solved by imposing the matching conditions v + (9f) = V-(9f) and v + (9b) = V-[2i r + 6b) where 9f{r) and 
6b(r) are the positions of the interface from quiescent to excited and back to quiescent, as a function of the radius r 
(see Fig. 1). Defining u> = co/g, g± = g±/g, the solution can be written as 


v±(r,Q) = g± + A ± [r)e* 


( 6 ) 


where 


A -(r) 


_ ( e -0//£> _ e -<V“) 

(#+ — 9- ) ^ ( 


A+(r) = A-(r) - ( g+ - g-)e 9f/ “ 


This can then be substituted into the “eikonal” equation (2) to find the interfaces d/(r), 6b(r). Near the tip (9 = 0, r = 
rt), smoothness requires —9b{r) = 9f(r) = a^/r — rt, for some positive a. This condition allows for the determination 
of the rotation frequency to. It is easy to check that the field value v takes the value </_ at the tip and for all r < rt. 
This solution has several unattractive features. Since the interface positions 9f b vary as y/r — r t , 


dv 1 

dr y/r - rt 


( 7 ) 


for all 9 when r is close to but above the tip radius r t . Also, ^ has a finite jump discontinuity at r = r t , 9 = 0, so 
that the derivative of v along the interface has a jump discontinuity at the tip where the front and back meet; this 
leads via the relation (2) to a similar jump for the derivative of the curvature. In addition, the normal derivative 
of v has a finite jump discontinuity across the interface. These difficulties have had the practical effect of making it 
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impossible to do a full stability analysis of the spiral solution jb|] and in general raises questions concerning whether 
the introduction of finite diffusivity dramatically alters the conclusions one derives via this construction. We now 
show that in fact finite diffusivity smoothes out this singularities without making any quantitatively significant change 
in the interface shape and selected rotation frequency. 

To proceed, we use an integral equation formulation of the problem ||. The field v — g± obeys a homogeneous, 
linear field equation with a fixed discontinuity across the interface. Following standard manipulations, we can therefore 
express v in the rotating frame as 


v±(r, 0) = g± + (<?+ - g-)D J ds' (h 1 • V'G - ^-h'gG^j 


where G is the Green’s function for the equation 

(DV 2 + w ^- M ) GM;r',0') = 
An explicit representation of G which will be used below is 


-6(r - r')5{6 - 6') 


( 8 ) 

(9) 


G(r, 9;r',9') 



dt' — (r ,2 +r' 2 

- t-) eX[>{ - 


2rr' cos {6 — O’ + u>(t — t'))) 
4 D(t-t') 


-n(t-t')} 


( 10 ) 


The integral in (8) is over the entire interface with arclength variable s' and n’ is the unit normal which points 
outward from the excited region; h'g is the component of n in the direction of the 9' vector. For a given interface, 
this construction gives a field which obeys the field equation and the continuity condition; a full solution can then be 
found in principle by substituting the field value at the interface into the eikonal condition and iterating the interface 
shape until the equation is satisfied. What we will do instead is evaluate this integral for the Pelce-Sun (PS) interface 
obtained at D = 0. We will show that this field is completely free from singularities but differs from the Pelce-Sun 
field solution only by terms which are small in the small D limit. 

Let us first focus on some point away from either the interface or r = rt . In this region, we will show how our 
formulation reproduces the known results of Pelce and Sun. Since G will never be singular, the n' ■ V'G in (6) is 
order D and can be neglected. The integral fi' g ds' can be replaced by ±dr' (for front and back respectively). Now, 
the integrals over r' (in eq. 6) and t' (in the definition of G) are dominated by the saddle point contributions coming 
from 


r' = r + 0 (VD) 


t' n = t + (6 — 6' — 2Trn)/to + 0 (y/D) 


If 9 > O ', n runs from 1 to oo; otherwise n goes from 0 to oo. Doing the Gaussian integrals around these saddle points 
leads after some algebra to 

v±=g±- (g+ - g~) J^exp(0 - Of - 2?vn)/u - ^exp(6> -9 b - 2irn)/ti (11) 

\ n n / 

where the n values in each sum obey the aforementioned rule. In the excited region, for example, Of > 0 > 0 b , yielding 


v+(0) =g+- (g+-g-)e* 


.-Of /a 


> ^ ^ ^ — 2'Kn/u 
n —0 


— e 


— 6 b /u 


1 ^ g-2* -n/ui 


( 12 ) 


with a similar expression obtainable for V-{6). Doing the sums, we directly recover the result of eq. (6) above. 

We now wish to understand the departure from the Pelce-Sun solution due to finite diffusivity. As we have seen, 
there are a number of different regions where the Pelce-Sun solution breaks downs and exhibits singular behavior. 
These breakdowns are all cured by boundary-layers whose width vanishes with D. Interestingly enough, as we shall 
see, the scaling of the boundary-layer width differs in each region. We first examine the region r ~ rt, but still away 
from the interface, i.e. 0 is not close to 0. We saw in eq. (8) that dv/dr was singular in this region. Examing our 
integral representation for v was see that the saddle point contribution from the t' integration is unchanged; however, 
the Gaussian integral with r' ~ r + 0 {V~D) cannot be blithely extended to infinity since the lower limit of the r' 
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integration range is r t . Also, we must use the fact that the integral of the Green’s function along the front and back 
have opposite signs and cancel to lowest order for r' ~ r*. Those considerations lead to the expression 

v~{9) = g-- {g+- g~)ur t \/D J dr' r'; 0, 9') | e - =0 (Of(r') ~ 0 b (r')) (13) 

where 

r = rt + y/~Dr, r' = r t + V~Dr'. (14) 

Plugging in the expression of the Green’s function and performing the t' integral, we obtain 

ap 

v~(0) =g- + 2(5+ - g_)aP 3/4 — 

with 


m = £ 


dr' 


: exp 


— (r 2 + f' — 2ff / ) 
4(27rn — 0)/w 


+ (0 — 2n n)/dj 


Jo \J — 9)/u> 

F can be expressed in terms of the parabolic cylinder function for index v = —3/2: 

«-(0) = 9 - + £± ^ £)1/4 W (E n ^ 1/4 exp(5(5 n /2) 

’ ex P (i£r) d - 3/2 )) 


(15) 


(16) 


with 5 n = 2(27rn — 9)/uj. This structure represents a boundary layer of (minimum) width y For negative 

f, the parabolic cylinder function decays as a Gaussian and u_(0) approaches </_; for large positive f, D grows and 
cancels the exponential factor, leading to the expected behavior V-(9) — 5 _ ~ y/fD 1 / 4 ~ -^/r — rt. Note that the 
change in the field itself is negligible as it vanishes in the small D limit as Z? 1 / 4 , even though the derivative is of order 
D~ 4 J 4 and approaches infinity. 

So, we have shown how integrating over the Pelce-Sun interface with the finite D Green’s function leads to a 
regularization of the infinite slope discontinuity at r = rt, without modifying the field. In fact, this lack of field 
modification is true everywhere including the interface. This follows from the fact that we have already shown that 
our field construction agrees with the PS field everywhere that is a distance 0(1) away from the interface and the PS 
fields are continuous across the interface. Field derivatives near the interfaces, on the other hand, are significantly 
modified due to the presence of diffusion-induced interfacial boundary layers. We now turn to a discussion of the form 
of these layers. 

Let us consider first the interface away from the tip. Expressed in terms of our discussion so far, the only necessary 
modification is to the n = 0 saddle point in the t' integral. Specifically, if 9 — 9f(r) =9 is small (i.e. we are close to 
the front), we need to keep an extra term in the argument of the cosine appearing in the Green’s function 

9 - 9f(r') + cu(t — t') ~ 0 + ^p-(r — r') + w(t - t') (17) 

or 

The modified t' saddle point occurs at 

v(t-l') = -0- §/f(r-r') (18) 

(which must of course be positive). Doing the t' integral leaves us with an expression of the form 


dr’ 


0+BL(r-r>) 


: exp 


~(r - r') 


■M2 


{^\ d -t(r-r’) + 9\) 


+ (0 + -^ 7 (r-r'))/u, 


(19) 


where the integral ranges over those values of r , consistent with the aforementioned positivity constraint on t — t'. 
This expression implies there is a crossover in behavior from the previous off-interface structure when 9 is order D 
and that the range of relevant r' variation is also of order D. The overall structure of the v field can be shown to 
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( 20 ) 


v = v ps + Dv 



and v is constant on one side, exponentially decaying on the other side, of the interface. This form is necessary to 
cure the finite slope discontinuity across the PS interface. 

In fact, the width of this interfacial boundary layer is proportional to |^f |. Since this derivative diverges near the 
tip, the boundary layer is much wider than order D in this region. To analyze this region, we modify Eq. © above, 
making explicit the square root dependence of Of with r' near the tip: 

9 — 9f{r') + ui(t — t') + u(t — t') (21) 

(where 6 is measured from zero angle). Letting r = r t + r, r' = r t + r' , we obtain for the G term 


-d 

~d9 


dr' 2a\fp 


47 tD(— 9 + ay/W) 


: exp 


— (f - f')- 


4 Dt-d + ay/fi) 


+ (9 — aVp)/u> 


( 22 ) 


The crossover from the previous r ~ r t behavior at finite 9 occurs at 

9 ~ V? ~ £> 1/3 


so that the boundary layer size in f is f' ~ D 2 ! 3 4 5 6 7 . 

Hence the leading boundary layer structure near the tip is 111 


= g-+DV 3 i 


f 


r — r t 


9 


\ £>2/3 ’ £)l/3 


(23) 


(24) 


Note that this is exactly the type of relationship that we need. The tangential derivative discontinuity of vps can be 
explicitly canceled by the 9 derivative of the second term; if we iterate the equation perturbatively, we will find that 
the curvature of the PS solution needs to be corrected by an amount of the form D 1/,3 /(s/D 1 / 3 ) for small arclengths 
and this boundary layer will compensate for the jump in ^ across the tip. On the other hand, the normal derivative 
\ rt is actually of order D~ x ! 3 and hence appears to diverge in the Pelce-Sun solution. 

To summarize, we have shown how to construct a singularity-free spiral field, starting from the Pelce-Sun solution, 
by including the effects of small diffusion constant. The fact that this construction goes through without difficulty 
proves that the Pelce-Sun solution gives the correct interface shape and concomitant rotation frequency in the £> —> 0 
limit. However, the demonstration that the normal derivative of the field near the tip is actually divergent means 
that the stability problem is subtle; since a perturbation to the interface will in general move the tip in the radial 
direction, this induces large field changes which must be explicitly balanced by interfacial structure on the boundary 
layer length scale D 1 / 3 . Work on this stability problem is in progress. 
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